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The fluxes of light antinuclei A< 4 induced near earth by Cosmic Ray interactions with the inter¬ 
stellar matter in the Galaxy and with the Earth atmosphere are calculated in a phenomenological 
framework. The hadronic production cross section for antinucleons is based on a recent parametriza- 
tion of a wide set of accelerator data. The production of light nuclei is calculated using coalescence 
models. For the standard coalescence model, the coalescence radius is fitted to the available experi¬ 
mental data. The non annihilating inelastic scattering process for the antideuterons is discussed and 
taken into account for the first time via a more realistic procedure than used so far for antiprotons. 
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I. INTRODUCTION 

Antiprotons in Cosmic Rays (CR) have been exten¬ 
sively studied both experimentally and theoretically over 
the last few decades, with the general purpose of mea¬ 
suring their flux and understanding their origin. It is 
now generally agreed that the dominant part of the 
p CR spectrum is a secondary flux originating from the 
hadronic production induced by CRs on the interstel¬ 
lar (IS) medium (ISM). This agreement is grounded on 
the ability of the calculations based on this assumption 
to reproduce the data. Recently some new prospects 
have been outlined, strengthening the motivations for the 
study of the p flux []J , and extending the interest for CR 
antimatter to other light antinuclei, antideuterons d in 
particular, with the emergence of new astrophysical is¬ 
sues. The p and d production in neutralino annihilation 
has been considered as a possible signature for the Dark 
Matter constituents in the universe [2|. Antiprotons and 
antideuterons have also been considered as evaporation 
products of primordial black holes (PBH) and their flux 
at earth calculated in the perspective of searching for a 
possible signature of the source in the experimental spec¬ 
tra IS A common feature of these studies is that the 
calculated fluxes are extremely small in intensity with 
their spectra peaked at low momenta. This requires the 
secondary galactic contribution to be accurately known 
for a significant search of the former to be considered. 

Another motivation for a careful examination of the 
antimatter production in the Galaxy is provided by the 
need of evaluating the galactic flux which will constitute 
a physical background for the forthcoming new experi¬ 
ments to search for primordial antimatter in the universe 
IIIS The evaluation can be performed using the 
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CR flux data from the latest generation of experiments 

Ilia 

In the atmosphere, The CR induced p flux has been 
measured recently [ll| at mountain altitude, and some 
antinuclei production is naturally expected to originate 
from the same source. These atmospheric secondaries 
will generate another physical background in the future 
experiments. Although it could be separated in prin¬ 
ciple from the galactic flux on dynamic and kinematic 
grounds, a soundly based knowledge of this background 
is required. Recently the phenomenology of the particle 
production induced by CRs in the atmosphere has been 
thoroughly investigated with the purpose of accounting 
for the large amount of new data available from recent 
ballon and satellite measurements mm- For the p flux, 
it has been shown in these works that this flux can be ac¬ 
counted for in a conceptually simple, albeit numerically 
elaborated, framework [lj, [T 3 ] . These calculations are 
extended here to the case of A>1 antinuclei. 

The primary purpose of the present study was to calcu¬ 
late the flux of antimatter particles of galactic and atmo¬ 
spheric origin, anchoring them as firmly and as widely 
as possible in the existing body of experimental data. 
The galactic and atmospheric production of p and d , 
as well as t , 3 He , and 4 He , and their flux in orbit 
near earth can be evaluated at so far unmatched level 
of confidence by using the available hadronic production 
data for these particles on the one hand, and the proven 
coalescence model on the other hand. A few prelimi¬ 
nary steps of this work have been covered already and 
published: The systematics of the p data over a wide 
range of incident momentum has been recently gathered 
and reanalyzed El , providing an accurate parametriza- 
tion of the inclusive antiproton production cross section. 
The first calculations of the atmospheric p flux were re¬ 
ported in [kJ- Other steps dealing with the coalescence 
models are quoted below. In addition, some important 
aspects of the production mechanisms of light composite 
nuclei and antinuclei relevant to the study are available 
from high energy heavy ion physics data, strengthening 
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its grounds : a key point is that the deuteron coales¬ 
cence parameter derived from empirical fits to the data 
has been found to be energy independent to within a few 
tens of percents, through a wide incident energy range 
including the range of interest here, for pA systems (see 
for example, and below). This feature is ex¬ 
pected to be valid for antideuteron coalescence as well, 
on general theoretical grounds. It is also supported by 
experimental indications El III 03 


The rescattering on the ISM (or on the atmospheric nu¬ 
clei) of the secondary particles in the propagation, and 
the contribution of this process to the low momentum 
flux of CR particles is an important aspect of the calcu¬ 
lations since several fluxes of primary origin and of major 
astrophysical interest are expected to contribute to the 
low momentum range, where the knowledge of the sec¬ 
ondary flux has then to be most accurate. The issue has 
been addressed in this work, with the same concern of 
keeping as close as possible to the experimental facts as 
well as following general principles. It is discussed in a 
dedicated section and appendix below. The subject of 
d rescattering has been addressed previously in ref [ 2 Cj • 


Finally, whereas the atmospheric antinuclei fluxes de¬ 
pend mostly on the proton and helium fluxes at the top of 
atmosphere (TOA), which are well measured, the galactic 
secondary antinuclei fluxes must be evaluated by means 
of a propagation model. Two classes of such models are 
widely used. The simplest one is the popular leaky box 
model (LBM), whereas a more realistic approach is the 
two-zone diffusion model (DM) [jj]. In the former, all 
the details about the transport conditions are discarded. 
These are globally accounted for by one single parame¬ 
ter, the confinement time of Cosmic Rays in the Galaxy, 
which, once fixed, leads to a balance between the produc¬ 
tion of secondary species (LiBeB, Z=21-23) from their re¬ 
spective primary parents (CNO, Fe) and the escape from 
the Galaxy, to reproduce the data [21| ■ This model gives 
correct results for the particle fluxes of stable nuclear 
species. 


The paper is organized as follows. The hadronic pro¬ 
duction cross sections and the coalescence model are dis¬ 
cussed in section im The total reaction cross sections 
used in the calculations for the antimatter absorption 
are presented in llll Al with the particles rescattering is¬ 
sue treated in section mm The antimatter propagation 
is discussed in section iTVl while the atmospheric produc¬ 
tion is presented in the following section m The results 
are discussed and the work is briefly concluded in the final 
section IvTH Some details on the physics of the rescatter¬ 
ing process are given and discussed in the appendix. 


II. ANTIMATTER NUCLEI PRODUCTION 
CROSS SECTIONS 

A. Production mechanism and coalescence process 

1. The standard coalescence model 


In the standard coalescence model, the invariant dif¬ 
ferential cross section for the (inclusive) production of 
composite fragments with nuclear mass number A, is re¬ 
lated to the inclusive production cross section of nucleons 
by a simple power law 123 : 
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where N a and N p are the A nucleus and nucleon pro¬ 
duction multiplicities respectively, and with p A = A ■ p* p . 
The coalescence coefficient B A can be expressed in terms 
of the coalescence momentum po and of the nuclear char¬ 
acteristics of the colliding system. 

In the case of deuteron or antideuteron coalescence out 
of matter-antimatter production, B A is defined by the 
relation giving the probability of finding a neutron and 
a proton within po distance in the momentum space l 22 l 
1231 , which leads to: 
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where 7 is the Lorentz factor. The multiplicity for parti¬ 
cle i is 


d 3 Ni 1 d 3 <Ji 
dp 3 a R dp 3 


where ctr is the total nucleon-nucleus reaction cross sec¬ 
tion. Assuming equal neutron n and proton p production 
cross sections and momentum distributions, the relation 
can be approximated as: 
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This relation can then be straigtforwardly extended to 
the coalescence of A nucleons as f23j : 
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Finally, the coalescence coefficient deduced from relations 
^ and ^ above, is given by: 


Ba = 

3 mt; 


( 5 ) 


The inputs for the calculation of the production cross 
section for mass A nuclei (antinuclei) then consist only of 
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the value of the Ba ( Ba ) parameter, the proton (antipro¬ 
ton) production differential cross section, and the total 
pp reaction cross section ctr. 

As quoted in the introduction, the deuteron coales¬ 
cence momentum or coalescence radius (see refs 0 ] for 
example for the definitions), derived from empirical fits 
to the data has been found to be practically energy in¬ 
dependent between a few hundred MeV per nucleon and 
158 GeV per nucleon incident kinetic energy (see e.g., 
mm and lammmmm for more detailed 
studies), as expected from the simple coalescence model 
f22l ]. The proper ty s eems to hold as well for mass 3 nu¬ 
clei and for 4 He j25[. The B ^-2 and 6^=3 coalescence 
parameters for p induced collisions on nuclei can thus be 
taken constant over this range. 

In the simple coalescence model successfully describ¬ 
ing the light nuclei production in pA and AA collisions 
at low and intermediate energies, the coalescence nu¬ 
cleons originate from the excited nuclear matter of the 
target (projectile) nucleus, and coalescence particles are 
produced in the target (projectile) rapidity range. In 
the case of high energy pp collisions, nucleon-antinucleon 
pairs originating from the hadronization process are pro¬ 
duced in the nucleon-nucleon (NN) center of mass (cm) 
rapidity region. The coalescence nuclei originate from 
these hadronization products which production cross sec¬ 
tion is maximum at rest in the cm and which therefore 
fly in the laboratory frame with the cm velocity. In these 
collisions, the rapidity distributions of the produced par¬ 
ticles is thu s sy mmetric with respect to the cm rapid¬ 
ity iilil In high energy A A collisions, 
both mechanisms are allowed, the rapidity distributions 
of light nuclei produced by nuclear coalescence of target 
nucleons are centered around the target rapidity, i.e., in 
the low momentum region in the laboratory frame for 
fixed target experiments, and symmetrically around the 
projectile rapidity (see [ll| for example). In these colli¬ 
sions, the coalescence production from single NN colli¬ 
sions is also expected to exist, with a much lower cross 
section however. 

There are therefore two sources of coalescence nucle¬ 
ons, produced in very different dynamical regimes: one is 
the production out of excited nuclear matter, the other 
is the production out of hadronization products from in¬ 
dividual NN collisions. The coalescence model which 
applies well to the first case doesn’t necessarily apply to 
the second. In this latter reaction however, the measured 
d and d production cross sections have been shown to be 
both consistent with a coalescence production mechanism 
as for nuclear systems [ 23 , 03 133, LMl H 3 , which indicates 
that the same approach can be used both for matter and 
antimatter production calculation. 

In addition, the similarity of the production mecha¬ 
nism for light nuclei and antinuclei was also suggested 
by various experimental indications in A A collisions. In 
this case however, the coalescen ce p arameter was found 
to be incident energy dependent fl5, 16 v 17, lSy38:). Note 
that the microscopic approach to the coalescence model 


was developed initially to explain this dependence of 
the coalescence parameter on the kinematic conditions of 
the reaction. For p+A systems however, the coalescence 
coefficients were found to be energy independent as men¬ 
tioned above, for d as well as for d 0, HU IMl ■ This will 
be confirmed below for d particles. 

These observed dynamical and kinematical features al¬ 
low therefore to make reliable predictions for the produc¬ 
tion in space of light antimatter nuclei from AA collisions 
(A standing here for any nucleus including proton), over a 
range of incident CR proton energies conveniently match¬ 
ing the useful range for the present study. 

2. Microscopic coalescence model 

A microscopic approach to the coalescence model has 
been developed recently in •3l| >4(}, which successfully de¬ 
scribed the mass 2 and 3 antinuclei production pA —>d X 
and pA —>3 X (3 standing for t and 3 He) experimental 
data. The model has been used in the present work as a 
complementary tool to the standard coalescence model. 

The two models have been used in the study, in order 
to provide an estimate of the theoretical uncertainty to be 
assigned to the calculations. This uncertainty was turned 
around however by the standard model being fitted to the 
available data. 


B. Antiproton production cross section 

The cross section data for the proton induced inclu¬ 
sive p production on nuclei and on the nucleon, have 
been reanalyzed recently, in the framework of a research 
program aiming at the understanding of the secondary 
particle flux produced in the atmosphere 0 - This was 
achieved using an improved analytical form used in pre¬ 
vious works for fitting the cross sections Cl , and it was 
shown that the data could be described with a fair ac¬ 
curacy over the incident energy range between 12 and 
400 GeV. The results of this work, i.e., for pA —> pX and 
pp —> pX cross sections, have been used in the present 
calculations, and its relevant part for the present study 
is repeated here for convenience. 

The parametrization of the inclusive pA —> pX cross 
section data is one of the two major ingredients of the 
present calculations, on which results and conclusions are 
built. The pA —*■ pX cross sections obtained from the 
analysis of 01 are estimated to be accurate to within 
about 20% above Tp «1 GeV. The low kinetic energy 
p data 0.5 GeV< T p <1 GeV are probably less accurate. 
They tend to be somewhat underestimated by the ana¬ 
lytical form used 0 . No fully reliable estimate of the 
low energy T p <0.4- 0.5 GeV accuracy can be assessed 
because of the total lack of data over this range and of 
the high sensitivity to nuclear medium effects. Symme¬ 
try considerations also constrain the low energy data (see 
the discussion in 0 ) for the pp system. Unfortunately, 
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these cannot be fully applied to the pA data because 
of the nuclear medium effects on the final state in the 
p production, which may more than significantly distort 
the spectrum. This issue would deserve a dedicated study 
by itself. For the pp —► pX cross section, two versions of 
the parametrization are available: the A=1 version of the 
pA parametrization (below referred to as I) , and a ded¬ 
icated simplified version of the formula with no nuclear 
mass dependence in the analytical form, which provided 
a markedly better fit to the p production in pp collisions 
data (below referred to as II). The two versions provide 
values of the integrated cross sections (multiplicity) dif¬ 
fering by about 20%, the [pA —> pX]a= i fit (I) providing 
the larger values. Both parametrizations have been used 
in the present work for the pp cross sections. 


C. Antideuteron production 

Galactic antideuterons are expected to be produced 
in pA and AA collisions dominantly by the elementary 
pN —► dX process via a coalescence mechanism out of 
the hadronic production of NN pairs as discussed previ¬ 
ously. Another competing production is possible however 
in space, from the pp —> dX reaction. This reaction has 
to be considered since although the p flux is much lower 
than the p flux, only one NN pair has to be produced 
in the kinematical domain where the ’’incident” p in the 
final state can coalesce with a n particle produced in the 
collision, providing a lower energy d flux than the other 
reaction but with a comparable magnitude. Both produc¬ 
tion channels are discussed in the next two subsections. 


1. Antideuteron production in the pA —* dX reaction 

The simple coalescence model described above has 
been applied to a set of antideuteron cross section pro¬ 
duction data available in the literature 13 HE 3 |36j. 
Fig. □ shows the results obtained using the coalescence 
formula Q] with the p spectra obtained as in 0 , and fit¬ 
ting the value of the coalescence momentum to a selected 
sample of data, with the result po=79 MeV/c. As seen 
on the figure the data are well reproduced by the calcu¬ 
lations over a range extending from 200 GeV/c incident 
momentum up to the ISR center of mass (cm) energy 
s 1 / 2 =53 GeV. In total 34 data points have been fitted. 
The results are given in table []] 

The 70 GeV data from refs [27, 3] have been discarded 
from the search. They were found not to be compatible 
with the other data, since they provided exceedingly large 
X 2 values. These data are listed in table HU It is seen 
in the table that the data measured on Pb target have 
a not so bad y 2 per point value (bottom line in the ta¬ 
ble). They have not been used in the analysis however, 
for the sake of consistency. Although the observed dis¬ 
agreement could be of physical origin since the discarded 
data have been measured at an incident energy below the 


momentum range of good agreement in Fig.^ they have 
not been treated on the same footing as the other data 
because the p yields measured in the same set of exper¬ 
iments were already found to be inconsistent with the 
measurements from other works 0- The comparison of 
these data with the standard coalescence model calcula¬ 
tions using a coalescence parameter fitting the data sets 
listed in table 0 are shown on Fig. |3 The disagreement 
clearly varies from moderate to sharp, with no clear trend 
pointing to some possible physical origin. The data are 
very overestimated or very underestimated by the cal¬ 
culations for the same incident energy for the different 
measurements. In the calculations of the d and 3 







FIG. 1: Comparison of the deuteron production cross section 
calculated with the standard (dashed lines) and microscopic 
(solid lines) coalescence models, with the experimental data 
from pA collisions at lab momenta 200 GeV/c on Be and Al 
targets (top left ) \?Al . at 300 GeV/c on Be, Ti and W targets 
(Top right) \2^L and at 53 GeV center of mass energy on the 
proton from Wa[ (full circles) and f 3A 1 (open circles, bottom 
left), and from \3f\l (bottom right). 


yields, the threshold effects were taken into account by 
means of the phase space of the final state as described 
in ref Eol. 


2. Antideuteron production in the pA —* dX reaction 

For the antideuteron production calculation in the 
pp(He) —> dX reaction(s) either in the simple or 
in the microscopic, coalescence model, the inclusive 
p production cross section from the same incoming chan¬ 
nel is needed. This cross section being unknown exper¬ 
imentally, an evaluation procedure has been used to es¬ 
timate it on the basis of reliable approximations. In 
the assumed production mechanism the pp —> dX re¬ 
action produces one single NN pair out of which the 
h (antineutron) will coalesce with the ’existing’ p in its 
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FIG. 2: Comparison of the antideuteron production cross sec¬ 
tion calculated with the standard coalescence model as de¬ 
scribed in the text, with the inclusive experimental data from 
pA collisions at 70 GeV/c from refs \ 27. .1 37, IT?|/. The upper 
left frame shows the data measured at a set of angles, CP (full 
circle, solid line), 12 mrad (inverted triangle, dashed line), 
27 mrad (full triangle, dotted line), and )7 mrad (open circle, 
dash-dotted line). See text for details. 


Experiment 

target 

Pine Of ^ 

GeV/c GeV 

N x 2 /N 

Bussiere et al. [2£] 

Be 

200 

5 

0.95 


Al 

200 

3 

1.7 

Cronin et al. [2T] 

Be 

300 

1 

4.4 


Ti 

300 

1 

1.3 


W 

300 

4 

13.7 

Alper et al.-, Gibson et al. [3h] 

P 

y/s = 53 

8 

1.3 

Albrow et al. [3b] 

P 

yfs= 53 

3 

1.0 

Armitage et al. [4£] 

P 

y/s = 53 

9 

4.6 


TABLE I: Experimental d production cross section data used 
in the coalescence model analysis to fit the coalescence mo¬ 
mentum parameter as discussed in the text. N is the number 
of experimental points used in the fit. 


final state, into a d , i.e., pp —> (pnjnpX —> dX. The 
following assumptions based on the knowledge of the 
pp —> pX and pp —> fiX processes, have been made for 
this evaluation: 

1. The probability to produce an antiparticle by 
hadronisation in a pp collision pp —> fiX in the 
final state energy range of interest, i.e., excluding 
the quasi elastic charge exchange processes to the 


Experiment 

target 

Pine 

(GeV/c) 

N x 2 /N 

Binon et al.; Antipov et al. [4F] 

Al 

70 

6 2.10 4 

Abramov et al. [27, jlT] 

P 

70 

2 1.10 3 


Be 

70 

3 3.10 2 


Pb 

70 

4 15.4 


TABLE II: Experiments whose results were not included in 
the search. The chi-squared per point values (rightmost col¬ 
umn) were obtained using the coalescence parameters fitted 
to the selected data listed in table Q] Same definitions as in 
previous table. 


resonance region, is the same as the probability to 
produce a p in a pp collision pp —> pX: 

En iz~(P p nX ) ~ E P~fr I ' ( dPP ~p x )- 

d p n d Pp 

For the pp(He) —> fiX inclusive cross section, pa- 
rameterisations II and I were used for pp and pHe 
collisions respectively. 

2. The pp —> pX cross section is assumed to be the 
same as the pp —> pX cross section, namely: 

„ d?(Jfi , ^ „ d 3 (Tf, , 

E Pj3^(PP pX ) ~ Ep ~dP^( pp pX ) 

This approximation makes sense since the collision 
dynamics are similar for the two systems at high 
enough collision energies. It is probably not correct 
however for low energy secondary ps from nuclear 
targets because of the dominance of the annihila¬ 
tion cross section in p N collisions at these energies 
and of the subsequent interactions of the produced 
ps with the nuclear medium. For the proton inclu¬ 
sive cross section, the parameterisation from li^ 
has been used. 

An additional correction has to be made for the en¬ 
ergy dependence of the d production cross section 
near the threshold, assumed to follow the phase 
space available to the particle |4£j . 

The differential cross sections for the pp —► dX and 
pp —> dX reactions calculated using the coalescence 
model(s), are compared on Fig. |3] for incident energies 
near their respective threshold. For the latter reaction, 
the calculation was performed within the approxima¬ 
tions described above. It can be seen on the figure that 
the pp —> dX cross section is much larger than for the 
pp —> dX reaction, the maximum values being about four 
orders of magnitude apart. The distributions are also 
centered at very different particle energies, correspond¬ 
ing to the respective NN center of mass velocities in both 
cases, and being thus much lower for the pp —> pp(nn) 
reaction than for pp —> pp(pp)(nn) reaction. These fea¬ 
tures clearly show that the pp — > dX contribution has to 
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FIG. 3: Spectral distribution of the antideuteron production 
cross section in the two reaction schemes pp —> dX (solid 
line) and pp —> dX (dashed line) considered, calculated with 
the standard coalescence model as described in text. Tp and 
T p are the production thresholds. 


be taken into account in the calculations of the galactic 
d flux. 


significant result for the present study which purpose is 
to fix the orders of magnitude of the corresponding sec¬ 
ondary CR flux. It is interesting to note on the figure, 
that the experimental production cross section is found 
significantly larger for t than for 3 He for all values of mo¬ 
menta where both have been measured. Should this dif¬ 
ference be assigned to Coulomb effects in the coalescence 
process is an open question which would deserve a ded¬ 
icated investigation. It must be noted however than the 
production cross sections for t and 3 He particles measured 
in the same experiments were found almost identical [28l . 




D. Mass 3 and 4 antinuclei production cross section 

The production cross sections of mass 3 antinuclei 
(t 3 He ) in pA collisions are extremely scarce. They 
have been measured in the past only at 70 GeV I44ll4a . 
and at 200 and 240 GeV i2Sj. 

In the Serpukhov experiments 0, the t and 
3 He production cross sections have been measured at 
small production angles on Be and Al targets. The p and 
d data from the same series of experiments have been 
discarded from the sets of data fitted in the present anal¬ 
ysis, as discussed previously. For the results of the mass 
3 antinuclei production cross section measured in these 
experiments, no consistent interpretation could be found 
either |4(|, albeit the standard coalescence calculations 
are found in reasonable agreement with the data, i.e., 
typically less than one order of magnitude larger. 

The production cross sections of mass 3 antinuclei have 
also been measured at CERN in a series of two experi¬ 
ments |2|j. The measured p and d cross sections have 
been discussed in the previous subsection. The mass 3 
data are compared on Fig. 0 to the calculated cross sec¬ 
tions in the coalescence models. The calculated cross 
sections are found within one order of magnitude from 
the data points for the two studied targets, with the two 
coalescence models providing values approximately fram¬ 
ing the data points, the microscopic model providing the 
upper limit. This agreement can be considered as fair as 
far as the order of magnitude is concerned. It is quite a 


FIG. 4: Comparison of the production cross section for mass 
3 antinuclei calculated by means of the standard (dashed line) 
and microscopic (solid line) coalescence model as described in 
the text. 


III. OTHER CROSS SECTIONS REQUIRED 
FOR THE CALCULATIONS 

A. Total antimatter reaction cross sections on 
nuclei 

The total reaction cross sections UR^aA) for light 
antinuclei a collisions on nuclei A are needed for compu¬ 
tation of the absorption term in the LBM or DM trans¬ 
port equations. 

1. Antiprotons 

For the galactic propagation, the (Jr(pp) total reac¬ 
tion cross sections for the pp collisions were calculated 
following ^8(, while for p 4 He collisions they were calcu¬ 
lated from the former in the Glauber approximation as 
described in the next section (see ^(J for more details). 

For the atmospheric propagation of ps in the atmo¬ 
sphere, the total reaction <jr and annihilation a ann cross 
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FIG. 5: Total p reaction cross section on the on 12 C (full 
triangles, see \If] for the sources of data), compared with the 
fit using relation(dashed line). The pC annihilation data 
(full squares) are compared with a fit using relational (solid 
line). The pp annihilation data (full circles) are compared 
with the fit to the data points from El/ (dot-dashed curve). 


sections for pA collisions on light nuclei ( 14 N and le O) 
were calculated using the following formulas (in mb): 


cr R = (257.8 + 


^) 4 )2/3 ( 6 ) 


a ann = 0.661(1 + 0.0036T -0 ’ 774 - 0.902T ool51 )A 2/3 (7) 

T being the particle kinetic e nerg y. This relation is based 
on a similar formula given in [43 (p 39 of this ref.), with 
the A 2 / 3 dependence added. The resulting fit is shown 
on Fig.0 


2. Antideuterons 

A few experimental data for the total reaction cross 
section cr R (d A) are available 0]. The same func¬ 
tional dependence on the nuclear mass number A as ob¬ 
tained in 491 has been used in this work: 


<jr = 105A 2 / 3 mb 


( 8 ) 


The total reaction cross sections have also been calcu¬ 
lated using the Glauber approximation and the parame¬ 
ters of the elementary pp scattering amplitude using the 
NN amplitudes from the analysis of [5lj]. The results 
are compared on Fig. 0 with the experimental results on 
the proton target and with the empirical parametrization 
from relation 0 See ^(| for more details. 


3. Mass 3 and 4 antinuclei 

For the mass 3 and 4 antinuclei, the total reaction cross 
section has been calculated in the same eikonal approach 
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FIG. 6: Top: Total dp cross section data from E3/ compared 
with the calculated values in the single scattering approxima¬ 
tion of the Glauber model, as explained in the text. Note the 
good agreement between calculations and data. Bottom: Com¬ 
parison of the total d 4 He reaction cross section obtained us¬ 
ing the empirical relation 0 (full circle, the error bar is from 
UAl ). with the Glauber approximation results (dashed line). 
The calculated total (full line) and elastic scattering (dash- 
dotted line) cross sections are also given, for completeness. 
See text for details. 


as for the p A and d A systems, using the known mat¬ 
ter distribution of the colli ding nuclei and the nucleon- 
antinucleon amplitudes |4(J |5ll ■ 


B. Non annihilating inelastic rescattering 

Before proceeding to the calculations of the propa¬ 
gated fluxes, the important issue of the rescattering of 
the transported particles is discussed in this section and 
more extensively in the appendix. 

The low kinetic energy range - T <500 MeV per nu¬ 
cleon - antideuteron flux is most sensitive to the possi- 




















ble primary d flux originating from exotic astrophysical 
sources such as primordial black holes @ or dark mat¬ 
ter (neutralinos) annihilation Q. Any calculations of the 
secondary flux aiming at a good accuracy for low energy 
particles must then take into account all the significant 
effects contributing to populate this energy-momentum 
region. Note that this argument holds qualitatively as 
well for antiprotons or any other particle or antiparticle 
propagating in the interstellar medium. The rescatter¬ 
ing of particles involves some energy loss and thus the 
transfer of a fraction of the flux from a given energy to 
a lower energy. The energy loss induced by the elas¬ 
tic scattering of the transported particles at the energies 
considered here involves only small momentum transfers. 
It is negligible. The inelastic process dp -^>dX may in¬ 
volve large energy-momentum loss of the scattered par¬ 
ticles. A reliable description of the secondary flux in 
the low energy region should thus take into account the 
component induced by the dp —>dX Non Annihilating 
inelastic (Re) scattering (NAR) reaction of the particles 
propagated in the ISM. This rescattering component will 
be referred to as tertiary in the following. 


1. Principles 

It has been argued recently that d (or symmetrically 
d) particles incident on a nucleon or on a nuclear target 
should have a small inelastic scattering cross section be¬ 
cause of their natural ’’fragility” originating in the small 
d (d) nuclear binding energy, since a deuteron bound by 
only 2.2 MeV easily dissociates in a collision. This intu¬ 
itive argument is misleading however, and appears not to 
be correct after a careful examination of the problem. It 
can be invalidated both on empirical evidence and on for¬ 
mal grounds. The former consists of the existing experi¬ 
mental evidence of deuteron [ 0 , [E 3 - and more generally 
nuclei 01 - induced nucleon excitation cross section (see 
also section IT 11 11 21 below!. The latter requires a few de¬ 
velopments given in the appendix where it is shown on 
general grounds that the inelastic scattering cross section 
of light nuclear systems is in fact expected larger than the 
corresponding (p,p') or (p,p') cross section, basically on 
the simple argument that it is driven by the NN(NN ) 
elementary cross section folded with the matter distribu¬ 
tion of the colliding systems. This is in agreement with 
the available experimental facts. 

Effects of the isospin selection rules 
This picture is somewhat blurred however in the case of 
the inelastic (d,d’) reaction on the nucleon since the in¬ 
elasticity is strongly inhibited at low incident energies 
by the isospin selection rule, the isovector excitations in 
the target nucleon being forbidden for this reaction. The 
deuteron is an isoscalar particle with isospin 0 which can¬ 
not induce isovector (i.e.,with isospin quantum number 
T=l) excitations in a nucleon (A resonances), in single 
step transitions. This has a strong inhibition effect on 
the total inelasticity of the reaction since the first ex¬ 


cited state of the nucleon is the (isovector) P 33 A reso¬ 
nance which in addition largely dominates the excited 
nucleon spectrum when the transition is permitted in 
inelastic processes like (e,e’) or (p,p’). No direct nu¬ 
cleon excitation is thus permitted in (d,d’) below the first 
(isoscalar) N* resonance in hadron collisions, namely be¬ 
low the 1.4 GeV (Roper) and 1.52 GeV 0 resonances 
in the nucleon. 

The overall inelasticity is not totally hindered how¬ 
ever since two step excitations via the nucleon (A) reso¬ 
nance^) in the deuteron are allowed and have been ob¬ 
served experimentally It3 . Ibfiil . 

2. Practical method 

In account of the complex interplay between the un¬ 
derlying NN cross sections and selection rules which gov¬ 
ern the dp inelastic cross section as discussed above, the 
practical approach used here has been based on empiri¬ 
cal grounds. The energy integrated dp —> dX cross sec¬ 
tion has been inferred from the experimental values of 
the cross section for symmetric system p d —► Xd which 
should be identical at the same center of mass energy 0 ] • 
Fig. m shows the experimental values for p d —> ( mr)p d 
(n multiplicity of produced pions) reactions as a func¬ 
tion of the incident energy 0] . The total inelastic NAR 
cross section ai n (pd) has been obtained by summing up 
the p d ( mr)p d cross sections experimentally avail¬ 
able, leading to (Tinipd) « 4 mb. No attempt was made 
to evaluate the (expectedly small) contributions of the 
channels not known experimentally. The overall evalua¬ 
tion is thus quite conservative. Note that these data also 
support the arguments in favor of a non negligible d ( d ) 
inelastic cross section, given in the previous section. For 



P„ (GeWc) 

FIG. 7: Partial p d —> ( niv)p d cross sections from E2/ used 
to evaluate the total inelastic NAR dp^dX cross section. 

the momentum spectrum of inelastically scattered parti¬ 
cles, the same functional form as in [0 was used to fit 
the data measured in these works, where the pp —> pX in¬ 
clusive cross section was found to be independent of the 
longitudinal momentum in the center of mass p*, and 
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could be described in the laboratory as: 


d 2 a(pp —y pX ) 

dpdVl 


Pp ^(E — (3p cos 8) 
2-rtpt E 


610 p 2 e 


Pt 

0.166 


(9) 

where 7 et /3 are the usual Lorentz factor and particle 
velocity and pt the transverse momentum of the particle. 
The integrated cross section was normalized to the value 
determined as above for a given d incident energy. 

The ansatz used in previous works for the energy dis¬ 
tribution of the secondary particles was based on the lim¬ 
iting fragmentation hypothesis |^3|. The form used for 
the energy dependence of the differential cross section 
da(E,T 0 )/dE ~ 1/T 0 just corresponds to a constant dif¬ 
ferential cross section over the energy range of the pro¬ 
duced particles between 0 and the incident kinetic energy 
To in the laboratory frame. Experimentally, the pp —► pX 
differential cross section was shown to be largely indepen¬ 
dent of the longitudinal momentum p) of t he p roduced 
particles in the center of mass Q (see also 01). Fig. 0 
shows the inelastic scattering spectra in the laboratory 
for 2 GeV and 3 GeV incident kinetic energy protons 
obtained in the approximation of jj59f (flat spectra) and 
used in this work, after 0|. The effects of the observed 
differences on the galactic flux of particles are discussed 
in the section HV B 21 below. 



FIG. 8: Spectral distributions of inelastically scattered pro¬ 
tons (antiprotons) from a hydrogen target in the approach of 
(horizontal lines) and used in this work for evaluating the 
NAR contributions of the antinuclei flux, for 2 GeV (dashed 
lines) and 3 GeV (solid lines) incident kinetic energies. 


IV. GALACTIC ANTIMATTER PRODUCTION 
AND PROPAGATION 

In the previous sections, the hadronic and nuclear 
physics issues related to light antinuclei production and 


interactions have been reviewed and updated. The 
present section addresses now the flux propagation of 
the produced antimatter nuclei in the framework of the 
leaky box model. The model is simple, widely success¬ 
ful, and highly useful, with a sound predictive power for 
use with purposes such as the present one. It rests on 
a single effective phenomenological parameter - the es¬ 
cape length A esc [ 6 (| - incorporating the physics (diffu¬ 
sion and convection) of the transport process. The leaky 
box model avoids lengthy discussions about what should 
be the value of the galactic transport parameters in more 
realistic models. It is an economical way to obtain reli¬ 
able results. As such, it is very well suited to address the 
impact of the cross sections obtained here and in 0 on 
the p and d fluxes. 

The errors on the 3 He , 4 He and to a lesser extent, 
d production, are the dominant source of uncertainty. 
Even for ps , the inaccuracy on the production cross sec¬ 
tion has been estimated to be larger than that related to 
the propagation parameters in [I|. For heavier antinuclei 
A, the production rate is proportional to Pq 1 - An un¬ 
certainty of Apo on the coalescence parameter leading to 
an uncertainty Ad’j on the antideuteron flux translates 
into an uncertainty (A — l)A<f>j for these anti-nuclei. 

The following subsections compare in details the p and 
d fluxes with previous results for the various contribu¬ 
tions from nuclear reactions, diffusive reacceleration, and 
solar modulation. The results derived in the simple LBM 
are then compared to more realistic calculations per¬ 
formed in the DM framework. The section also includes 
the presentation of 3 He and 4 He fluxes and ends with a 
summary of the calculated fluxes. 


A. The Leaky Box Model for Cosmic Ray 
propagation 

In this model, the flux & A (T A ) °f the antinucleus A at 
kinetic energy T 4 is given by: 


$a(t a ) 


^esc (T A )X int (T A ) 

£?ISM [^esc (Ta) + ^int CT*)] 

x^[Q s a(Ta) + QT(t a )\ ( 10 ) 


with QJ c {T A ) and Qf(T A ) being the secondary and 
tertiary source terms discussed below, respectively, and 
where X- mt (Ti) = {m)/(a^ +1 ™(T A )) and X esc (T A ) are 
taken from |61|. Following |63], we use the quantities 
(m) = 2.05 x 10~ 24 g for the average mass of the in¬ 
terstellar gas, qism = 2.28 x 10 - 24 g cm -3 for the ISM 
average density and (<7£ +ISM (T J 4)) stands for the average 
reaction cross section on this gas. The ISM composition 
used is H : He : C : N : O = 1 : 0.1 : 5 x 10~ 4 : 
8 x 10 -4 : 8 x 10 -5 cm -3 . Actually, these numbers as 
well as the average density are not perfectly known. The 
requirement to fit the B/C ratio for a given choice for 
these quantities, leads to a peculiar parametrization for 
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Aesc(T^). However, a different choice than the above, if 
it would probably affect the normalisation of the fluxes, 
would have a minor effect on the shape of the distribu¬ 
tion. It is thus sufficient to check that the choice made 
for (m), Qism and A esc leads to the correct normalisa¬ 
tion of the spectra. Tested on the p spectrum, the same 
set of parameters can then be safely applied to the other 
antinuclei fluxes. 

Concerning the secondary Q s ^ c and tertiary Q)| r source 
terms in relation E3 above, the former results from the 
net creation of the antinucleus A from CRs interaction on 
the ISM (see Sec.[Hjl, whereas the latter results from the 
energy redistribution through inelastic non-annihilating 
(NAR) reactions of the produced antinuclei, discussed in 
details in Sec. mm above. Note that neither the ionisa¬ 
tion losses nor the reacceleration process were included in 
the LBM calculations. They were taken into account in 
the DM results presented in Sec. IIV B 3 cl however. The 
main new features of the present calculations with re¬ 
spect to previous works are: 

1. The anti-nuclei production cross sections used are 
more tightly constrained by the available data M- 

2. The standard coalescence model for the produc¬ 
tion cross section of composite particles has been 
calibrated and validated on the available inclusive 
production cross section of light antinuclei (d , t , 
3 He ), and an independent microscopic (diagram¬ 
matic) coalescence approach was also used. 

3. The inelastic rescattering of the transported parti¬ 
cles are taken into account in a more realistic way 
than in previous works. 

4. The contribution of the pp —> dX reaction cross 
section was included in the calculated d flux. 


B. Antiprotons and antideuterons 


For the p flux, the two source terms are 


p,He,CNO H,He,CNO 

Q sec (T p ) = 2 £ Y, 

i=CRs j=ISM 

r °° da i+ i 


xAirrin 


, r i r r_ v p 

' 6 m p u ± p 


(Tp, Ti)^i{Ti)dTi (11) 


Q tei (T p ) = 47T.7 


f°° drrPP^P x 

2 / (T' p ,T p )%(T^)dT' p 


It dT p 


-2<r PX (Tp)$p(T p )\ (12) 


where n p is the hydrogen number density in the ISM in 
cm -3 , da 1+ )/dT p is the differential p production cross 
section, daP p ^ x /dT p is the differential inelastic non¬ 
annihilating cross section for ps with incident energy T p 
emerging from the collision with an energy T p < T'-. and 


°"frfa^ 7)A ^e total inelastic scattering pp —► pX cross 
section. This latter term involves the quantity <P p {T p ) in 
the integrand, requiring a numerical method to be used 
to solve Eq. W- Following [3 , the tertiary contribution 
on He is taken into account assuming a mere scaling fac¬ 
tor 4 2 / 3 to the corresponding cross section on Hydrogen. 
The factor of 2 in front of both terms in 1 1 II takes into ac¬ 
count the n (decaying into p ) production cross section, 
assuming the p and n production cross sections to be the 
same. 

For antideuterons, the source terms are similar (with 
p labels changed to d) with only the pp , pHe and Hep 
incoming channels for the production reaction (thresh¬ 
old 16m p for the pp incoming channel) being taken 
into account (heavier components neglected). The new 
contribution Q tel (pp and pHe, threshold 6 m p ) which 
was assumed to be negligible in [ 2 ] was included here. 
In Q sec (T£j the term da CRs+1SM /dT^ was evaluated by 
means of the coalescence models discussed in Sec. rrm 

Note that for the primary fluxes of p and He, the 
parametrization provided by the AMS01 experiment was 
used. Since the flux measurements of the BESS [§] and 
AMS01 U experiments are now compatible, the primary 
fluxes can then be considered as being a minor source of 
uncertainty of the calculations Q . 


1. The secondary source term 

Fig. EH shows the various contributions to the source 
term Q sec (T p ) as defined in relation IT21 (independent of 
any propagation model) for the antiproton (left panel) 
and antideuteron (right panel) fluxes. 

For the antiproton flux, the pp reaction contributes to 
about 56% to the p production, the pHe to 24%, Hep 
to 12% and the He He reaction for 6%. The latter has 
been evaluated by scaling the pp —► pX cross section 
using the same procedure as described in .£3]. The re¬ 
actions p — CNO and CNO — p contribute to less than 
2%, while He — CNO and CNO — He have been ne¬ 
glected in account of the very low CNO nuclei flux and 
IS density. These components are compared on the fig¬ 
ure to those obtained in mu using the cross sections 
calculated using the DTUNUC event generator. At the 
production peak, for pp and pHe reactions, the calcu¬ 
lated yields are similar up to 10%. There is a trend of 
the present calculations to predict smaller cross sections 
than DTUNUC for particle energies above a few GeV. In 
the low energy range T^in <1 GeV, it is seen on the figure 
that DTUNUC clearly predicts a larger Jgsj or even much 
larger |l| target compositness effect than in the present 
work, with a larger low energy cross section calculated 
for composite targets (A>4) than for the proton target. 
This difference would deserve a further investigation both 
theoretical and experimental, of the low energy p yield 
in nuclear collisions. Note however that for the pp in¬ 
duced yield, the present calculations predict a larger flux 
at very low energies than the DTUNUC based calcula- 
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FIG. 9: Source terms for antiproton and antideuteron particles as defined in relation m Left panel: Q sec for the antiproton 
production from the present work (solid lines), compared with the results from two different papers using the DTUNUC event 
generator (jjj]: dashed lines and \hAl : dotted lines). Note that the curves taken from \6Al (dotted lines) correspond to a slightly 
different H and He input fluxes compared with the two other models presented. Right panel: Q sec for the antideuteron production 
for the different channels considered in this work. See text for details. 


tions. It must also be remembered that this energy range 
is marginally within the DTUNUC domain of validity. 

For the antideuterons (right panel), the contributions 
from the pp, pHe , and Hep collisions are peaked around 
4 GeV/n and the obtained distributions display a simi¬ 
lar bell shape as obtained for the IS p flux. This feature 
can be easily explained qualitatively since the calculated 
distributions of the d (p ) fluxes are driven on the low 
energy side by the rapidly raising d (p ) production cross 
section with energy above the d (p ) production thresh¬ 
old (see 1^3 f° r example), while the high energy decay of 
the distribution is determined by the rapidly decreasing 
incident CR proton flux with energy, folded with the nat¬ 
ural decrease of the high energy production cross section 
(see Fig®. 

These dynamical and kinematical characteristics (see 
snbsection lll C 21 are responsible for the pp —> dX contri¬ 
bution to the flux to be clearly dominant over the entire 
low energy range below 1 GeV/n with a value larger than 
the value reported in Sil by about one order of mag¬ 
nitude (note that the p 4 He —► dX contribution can be 
neglected here, as it can be seen by comparing curve A 
to curve B in Fig. ®. At energies above 2 GeV, this flux 
becomes rapidly negligible. 

Note that for this channel, the p flux was calculated 
conservatively using the parametrization providing the 
smaller values for the cross section. The range of in¬ 
tegration for the d flux calculation has been limited to 
100 GeV p incident energy (while for pp —> p + X, the 


integration was performed up to a few tens of TeV p 
energy to calculate the p production), since the p flux 
beyond 100 GeV is vanishingly small. The major source 
of error here is due to the uncertainty on p flux in the few 
GeV range. For example, a roughly twice larger d flux 
could be produced in this channel by taking the mean val¬ 
ues of the experimental data points for this flux, rather 
than the p flux obtained from the p data fits. For the 
other channels, the uncertainties are related to those for 
the p production, combined to the contribution induced 
by the coalescence model. It must be emphasized how¬ 
ever that, although the procedure used to evaluate this 
cross section is based on the underlying physics, the er¬ 
ror induced by the approximations made in this approach 
cannot be accurately evaluated. 

In conclusion of this section, it has been shown that the 
secondary pp —> dX term does contribute significantly to 
the IS galactic d flux, and that it is even dominant in the 
low kinetic energy Tj < 1.5 GeV/nucleon range. 

2. The tertiary source term 

The only tertiaries considered here are those created 
by the non-annihilating inelastic scattering of secondaries 
on the ISM (NARprocess).This contribution was first 
included for ps in [65( ■ 

The sharp kink observed in the low energy d spectrum 
in Fig. um right, can be better understood by examining 
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FIG. 10: Comparison of tertiary contribution using the standard 1/T parametrization (A) and the parametrization used in this 
work (B) for the p (left) and d (right) fluxes. The B tertiary component is further decomposed into individual contributions 
from three energy bands (labels 1, 2 and 3). See text for details. 


how inelastically scattered antiprotons or antideuterons 
are redistributed, since the redistribution mechanism 
may have significantly different effects on the respective 
shapes and intensities for the p and d flux. Since the 
Q ter term in Eq. 1121) has to be integrated numerically, 
the same iterative procedure as proposed in [I] was used. 
For antiprotons, a few iterations are required, while only 
one single iteration is needed for antideuterons. This is 
simply due to the much larger relative inelastic (NAR) 
cross section compared to total cross section for pp than 
for dp collisions (see section UlI Bl and appendix), which 
makes the numerical convergence slower in the former 
case. To better understand the details of the transfer 
process to lower energy, the tertiary contributions from 
three kinetic energy bands (T<2 GeV, 2<T<10 GeV, 
10<T<100 GeV) are displayed separately on Fig. EH 
The corresponding curves were obtained by first evaluat¬ 
ing the equilibrium spectrum, and then computing from 
this spectrum the tertiary yield from the chosen energy 
bands. 

For the antiprotons (Fig.^J left panel), the NAR cross 
section is large and the low energy tail is largely replen¬ 
ished. Several iterations are required to obtain the equi¬ 
librium flux allowing the second order NAR contribution 
to be significant, the first iteration replenishing both the 
medium and low energy bands, the former of those re¬ 
plenishing the low energies in the next iteration step. 

For the antideuterons, the NAR cross section is much 
smaller than for ps. The first consequence is that the 
tertiary component is almost two orders of magnitude 
smaller than the peak value of the secondary component 
(while it is less than one order of magnitude smaller for 
ps). One single iteration is thus required in the process 
and the contribution of second order interactions is neg¬ 
ligible. This explains the sharp upturn in the d flux that 
is not seen in the p spectrum. The secondary flux drops 
rapidly at low energy and the tertiary and secondary 


components become comparable only at s=s 300 MeV/n. 
In this case the NAR process doesn’t accumulate parti¬ 
cles at low energies as efficiently as in the antiproton case 
because of the much smaller NAR cross section (compare 
the two curves labelled 3 in the left and right panels). 

The effects of the NAR spectrum on the p and 
d distributions are also illustrated on the figure, with 
the curves labelled A and B corresponding to the 
parametrization used in |5fj and here (see section lTlI B 21 , 
respectively. Although for ds the results are hardly dif¬ 
ferent at the lowest energies, for ps , the more realistic 
inelastic spectrum leads to a low energy NAR flux larger 
by almost a factor of 2 than the other option. 


3. Discussion of the interstellar fluxes 

a. Comparison to previous published fluxes 

Left Fig.llllshows the results from the present work for 
the unmodulated antiproton fluxes compared with some 
previously published results. For the present calculations 
(full lines), the two curves correspond to the two different 
parametrizations used for the p production reaction in pp 
collisions M ■ The present results are close to those from 
|66j (dashed line). Both rely on better fits to the data 
than in ||3] ■ The latter, used in the present approach, 
would give a flux standing roughly midway between the 
two full lines for p energies above 1 GeV. The so defined 
range can thus be considered as the uncertainty due to 
the production cross section. 

In the results of the diffusion model :l), the flux is 
somewhat higher at high energy, while at low energy the 
effect of reacceleration is responsible for the much larger 
predicted flux (note the enlarged vertical scale in Fig. IIII) . 
However, it will be seen below that the two results are not 
so different after modulation. Note that as the parame- 
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FIG. 11: Partial calulations of Galactic antiproton (left panel) and antideuteron (right panel) interstellar fluxes for comparison 
purpose. Left panel : The present calculations (solid lines) are compared with those from Q/ (diffusion model with reacceleration, 
dashed lines) and from w6 J (leaky box model, dash-dotted line). The two curves from the present work (using standard NAR), 
obtained for the two different parametrizations of the production cross section, as well as those from fl]/, indicate the level 
of uncertainty of the results due to hadronic cross sections. Right panel (no tertiaries): Present results (pp, pHe and Hep 
contributions only) using the standard coalescence model (fitted to the d data) with po = 79 MeV (dash-dotted line) and 
the diagrammatic approach of coalescence (solid curve). The two lower curves show the present calculations using the same 
coalescence parameter po = 58 MeV (dashed line) as in 0/ (diffusion model, no reacceleration, dotted line). See text for details. 



JSnebc energy (G«V) 



FIG. 12: Left panel: IS p flux including all contributions, the two lines showing the uncertainty due to nuclear cross sections 
(see text). Right panel: d flux for the standard coalescence model (reference calculation, lower curve) and for the microscopic 
coalescence model. 


ters used for the present calculations provide the correct 
magnitude for the calculated fluxes (see left Fie . 11411 . they 
can be used confidently to evaluate the galactic fluxes for 
other antinuclei. 

Right Fig. |TT| shows the calculated antideuteron flux 
(dotted line) evaluated as in ref @] (dashed line), i.e., 
including neither the tertiary term (NAR) nor the pp in¬ 
cluded and using a coalescence parameter pg = 58 MeV. 
Note that the p and He galactic fluxes used in (H were 
slightly different however from the more recent measure¬ 


ments used here. Note also that in [gj, a diffusion model 
without reacceleration was used. The two calculations 
appear to be compatible. At low energy, the differ¬ 
ence observed in the antideuteron flux from the two ap¬ 
proaches originates from: 

1. The difference between the low energy p production 
cross section used here and that from [43 (domi¬ 
nant effect), and 

2. The kinematics of the p production near the thresh- 
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old, also modelled differently here than in i , pro¬ 
viding a larger antiproton - and thus antideuteron - 
cross section, over this range. 

The upper two curves (solid lines) on the figure, cor¬ 
respond to the standard model (reference calculation, 
lower curve) and to the microscopic model calculations 
discussed in Sec. HE! (upper). 

b. Full calculation 

Fig.ua displays the full calculation results for the IS 
fluxes obtained in this work. The left panel shows the two 
calculations showing the uncertainty on the pp produc¬ 
tion channel induced by the two parametrizations used 
(dedicated pp (II), and pA with A=1 (I), see section ITT bI 
and El). 

It was seen in sections fiTcl and in ITV R 1 1 that the un¬ 
certainty for the other channels that contribute to almost 
half of the flux, is similar, their contributions having been 
evaluated using widely proven calculation methods based 
on existing experimental data. Comparing with the re¬ 
sults obtained using the DTUNUC generator, large dif¬ 
ferences have been found for the various individual chan¬ 
nels contributions, although a rough agreement for the 
overall secondary production has been obtained. 

For the antideuteron flux (Fig. El right panel), the 
lower curve shows the reference results obtained with the 
standard coalescence model (based on a fit of d data) 
and including all the secondary and tertiary components 
discussed previously. The overall uncertainty associated 
to the reference calculation is of the same order of mag¬ 
nitude as for the p flux, plus the coalescence model con¬ 
tribution to be added (see section urn . It is estimated 
to be better than a factor of 2. The upper curve on the 
figure shows the microscopic coalescence model calcula¬ 
tion (which overestimates the experimental data on the 
average). 

At low energy, the pp —>c? X reaction can be estimated 
to be correct to within a factor of 2 regarding the accu¬ 
racy of the inclusive p available data. The tertiary term 
corresponds to a pretty conservative lower limit since, 
based on the evaluation method used (see section lHI D 21 . 
the upper limit being less than a factor of 2. From this 
point of view, some direct measurements of the pp —>d X 
production reaction, would be extremely useful to elim¬ 
inate the uncertainties corresponding to the approxima¬ 
tions made here. 

c. Comparison with the diffusion model 

Fig. ua shows the diffusion model calculations for 
the d flux, with the same common inputs as in the 
present calculations, i.e., samep and He galactic flux and 
hadronic cross sections, including the pp —>d X contribu¬ 
tion. The top panel shows the results without diffusive 
reacceleration for the DM calculations. The agreement 
between the two results is good, except in the 500 MeV/n 
range where the contribution of the pp —>d X contribu¬ 
tion is more salient in the present calculations than in 
the DM approach, making the predicted flux larger by 
a factor of about two. The middle panel shows the DM 
calculations including the diffusive reacceleration term 





FIG. 13: Comparison of the present IS (unmodulated) d flux 
calculations (solid line) with the diffusion model results 
(dashed line), without (Top) and with (Middle) reacceleration 
term included (see El) , and with modulation and reacceler¬ 
ation included (Bottom). The two lines for the DM results 
on the three panels correspond to the same estimated range 
of propagation uncertainty as evaluated in na and used for 
p in 0. 
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compared with the same present results as above. The 
observed effect of the reacceleration on the spectrum is 
similar to the effect of the solar modulation (compare 
with curve 1 in Fig. 1141 rightl. making the distribution 
at low energies below 1 GeV/n power-law shaped. The 
bottom panel shows the effects of the combined reaccel¬ 
eration term and solar modulation for the DM model, 
compared with the present calculation including the so¬ 
lar modulation effects. In this case, the general shape 
of the results are similar, the very low energy flux being 
predicted larger from the DM calculations than from the 
present LBM by a factor of 2 to 4. 


4- Solar modulation effects 

It has been seen in the previous sections that quite sig¬ 
nificant differences between the antimatter fluxes are pre¬ 
dicted in the astrophysically sensitive low energy range 
of the studied particles. In this section, the solar modu¬ 
lation effects on the calculated flux are incorporated and 
their influence on the various components of the calcu¬ 
lations together with the resulting overall accuracy are 
discussed. 

The situation for the low energy spectrum is drastically 
different for the modulated fluxes, i.e. at TOA, than it 
was before for the IS spectra. Fig. E] left shows the 
antiprotons flux using the standard, i.e. flat spectrum, 
NAR cross section (solid lines) and the parametriza- 
tion used here (dashed lines) for a solar modula¬ 
tion parameter (j> = 500 MV. Whereas, these two calcu¬ 
lations produced very different IS fluxes (see left panel on 
Fig. eg. they are hardly different once modulated. This 
is not surprising since the effects of the solar modulation 
are well known to occult the low energy range of the IS 
spectrum. This effect results in the present calculation in 
improving the agreement (up to a sound normalization 
factor however) with the results of the diffusion model 
with reacceleration (dotted lines) from [IJ. For each case 
the two bands correspond as before to the uncertainties 
associated to the hadronic cross sections involved: for 
the present calculations (solid lines) they include a) For 
the pp incoming channel, the uncertainties related to the 
I and II parametrizations; b) For pA (A >1) systems, 
the uncertainties of the fit to the data as quoted above 
(see J3 )i c ) For AA (A >1) systems those related to 
the (Glauber) calculation procedure and its ingredients 
EfUl . In [lj the uncertainty band was evaluated from the 
DTUNUC uncertainties with all the contributions of the 
relevant A A collision systems considered, combined to¬ 
gether. 

The d flux naturally suffers the same source of un¬ 
certainty due to the hadronic cross sections as for the 
p flux. Fig. right illustrates the status of the var¬ 
ious components of the d flux. The three curves la¬ 
belled 1, 2 and 3, show the relative contribution of the 
various source terms for the same modulation parame¬ 
ter cf) = 500 MV. As it could be expected from Fig. 0 


the p induced d production ( Qp ism ) dominates over 
the d rescattering (NAR) contribution ( Q ter ) at ener¬ 
gies smaller than about 500 MeV/n, where the latter is 
clearly negligible. This is even more true at higher modu¬ 
lation level (compare solid line, dotted line 1, and dashed 
line). 

Finally, the thick solid horizontal line on the figure 
shows a lower limit for the estimates of primary source 
- supersymmetric particles annihilation (SUSY) or PBH 
evaporation - contributions which fluxes in both cases 
drop rapidly to negligibly small values above 1 GeV/n. 
Some of these predictions Si give fluxes that could 
be one order of magnitude larger than the quoted limit. 
Hence, this confirms that the antideuteron signal is a 
probably good tracer to look for SUSY, likely better than 
antiprotons, provided that ds can be discriminated from 
ps and electrons at the appropriate level of selectivity, 
i.e., with a rejection power of the order of 10 5 , in the 
forthcoming experiments, which is quite an instrumental 
challenge. 

5. Uncertainties on the calculated fluxes 

In summary, in the context of using the p flux to derive 
constraints on new (astro)physics (primary exotic com¬ 
ponent), it can be considered that the NAR cross section 
is now evaluated with an improved accuracy, even if this 
accuracy is clearly not a major issue here since it has been 
seen that for rigidity R< 0.5 GV of the IS flux, all differ¬ 
ences are swept away by the solar modulation effects. On 
the other hand, it must be emphasized that the signifi¬ 
cant differences observed between the results obtained 
from data fits, and used here for the production cross 
sections, and those obtained from the DTUNUC predic¬ 
tions (for A A collision systems), would certainly deserve 
a more complete investigation sinceit is a major source 
of uncertainty for the flux calculations. 

Below 0.5-1 GeV/n, the contribution of the pp induced 
d flux which suffers many flaws as emphasized in II V B H 
is probably large. Within and below this region of energy, 
several sources of large uncertainties due to, in decreasing 
order of magnitude, hadronic cross sections, solar modu¬ 
lation, and propagation, have been seen to exist. These 
would dramatically increase the difficulty of reliable pre¬ 
cision astrophysics measurements over this energy range. 
Above 1 GeV/n, the dominant source of uncertainties 
remains the hadronic cross sections. 


C. Mass 3 and 4 antinuclei fluxes 

The same calculations have been performed for the 
propagated flux of mass 3 and 4 antinuclei using the same 
coalescence model, with p 0 = 78 MeV (note that the dia¬ 
grammatic approach used for the antideuteron case can¬ 
not be applied for A=4), and the same hadronic cross 
sections. Unlike antideuterons for which the pp reaction 
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FIG. 14: Modulated antiproton fluxes. Left panel: Present calculations, with the two linetypes corresponding to the single 
pp —* pX reaction, without (solid lines) and with (dashed lines) the NAR contribution. They are compared to the DTUNUC 
calculations (dotted lines), the two curves for the latter giving the overall production uncertainties as evaluated in )1] (see text 
for details). Note that the differences due to reacceleration, observed in the IS fluxes, are completely smoothened by the solar 
modulation effects. Right panel: Same calculations with the curves labelled 1-3 giving the fluxes obtained, for the full calculation 
(label 1), with the contribution of the new ingredients Q pISM switched off (label 2), and with both Q pISAI and the NAR Q ter 
contribution switched off (label 3), for the same modulation level (4>=500 MV). The (typical) level of primary fluxes evaluated 
by various authors is shown by the thick solid straight line on the left ElH/ (see text). 


provides a sizeable flux, similar reactions (e.g., dp —> dX ) 
are expected to lead to negligible contributions for heav¬ 
ier antinuclei production. This is because of the ex¬ 
pected smallness of the corresponding cross sections, due 
to the large energy-momentum transfer required to pro¬ 
duce antinucleon pairs which would breakup the incident 
(anti)deuterons. 

For the tertiary flux of mass 3 antinuclei, the NAR 
cross section would be large since the restriction to 
isoscalar transitions discussed for the d projectile would 
not apply (while it would for 4 He ), and a similar proce¬ 
dure as used before for ds could be applied. The exercise 
would be irrelevant however, in account of the smallness 
of the expected coalescence cross sections and of the cor¬ 
responding NAR flux. 

The uncertainties associated to these calculations are 
driven mainly by the uncertainties on the hadronic cross 
sections, boosted by the coalescence model exponentia¬ 
tion. For the 3 He flux, it is estimated that the calcula¬ 
tions are accurate within roughly one order of magnitude, 
which is a good enough level of accuracy for the purpose 
of this work. 


D. Summary for the calculated galactic flux 

The spectral distributions of the calculated galactic 
fluxes for mass 1 to 4 antinuclei are displayed on Fig. 1151 
(full calculation, standard coalescence model). The cal¬ 
culated d , 3 He 4 He flux have been multiplied by 10 4 , 
10 s , 10 12 respectively for presentation purpose. These 
fluxes are significantly higher than those derived in (20j, 


the difference being mainly due to the larger value of 
the coalescence momentum derived in the present work, 
the other smaller differences having been discussed in the 
text. The roughly twice larger d flux between the old and 
the new coalescence momentum translates into a factor 
4 for 3 He and 8 for 4 He . The upper of the two dashed 
lines for 3 He includes the addition of the t flux, this lat¬ 
ter nucleus decaying into 3 He with a half life of about 12 
years. 

V. ATMOSPHERIC PRODUCTION AND 
PROPAGATION 

A CR particle crosses a sizeable amount of the earth 
atmosphere before being detected by a balloon borne 
running experiment. The grammage seen by the CR 
in this process can be of the same order of magnitude 
as the grammage seen during its wandering through the 
Galaxy in a few tens of Myr (typically ~ 10 g/cm 2 at 
1 GeV/nuc). The same mechanism which leads to the 
secondary galactic antimatter flux, leads as well to a sec¬ 
ondary flux of the same particles via the interaction of 
CR particles with atmosphere (see 01 f° r a quantitative 
comparison). 

The CR induced p flux in the atmosphere has been 
measured recently at mountain level, where it is expected 
to be of purely atmospheric origin EH- Light antin¬ 
uclei have also been searched by the same experiment 
E. These atmospheric secondaries will constitute a back¬ 
ground for the CR flux of these particles in future em¬ 
barked experiments. For balloon borne experiments this 
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FIG. 15: Galactic flux for p (solid line), d (dashed), 
3 He (dash-dotted), and A He (dotted) antimatter particles. 
The lower (resp. upper) dashed line correponds to the case 
where the 3 H production is not taken (resp. taken) into ac¬ 
count (see text for details). 


background, produced in the atmosphere above the de¬ 
tector, must be evaluated to correct the measurements 
|68j . Besides, atmospheric particles can be trapped by 
the earth magnetic field and detected by satellite exper¬ 
iments. Although they could be separated in principle 
from the galactic flux on dynamic and kinematic grounds 
|69| . a basic theoretical knowledge of this background is 
required. 

As quoted in the introduction, the phenomenology of 
the particle production induced by CRs in the atmo¬ 
sphere has been thoroughly investigated recently with 
the purpose of accounting for the large amount of new 
data available from recent ballon and satellite measure¬ 
ments. The calculation of the antiproton flux was part of 
this effort [l3j. The latter have been recalculated in the 
present work with the tertiary (NAR) production taken 
into account, and extended to the case of A>1 antinu- 
clci production on the same dynamical basis as for the 
galactic flux, as discussed below. The present results 
complement and update our previous calculations 0] . 


A. Antimatter flux 

The atmospheric production of antimatter has been 
calculated along the same lines as in the previous calcu¬ 
lations of the atmospheric particle flux by the authors 
0 E7]. It consists of a 3D Monte-Carlo simulation 
program processing the propagation and interactions of 
charged particles in the earth environment 00- For 
antinuclei A < 4, the antimatter production cross sec¬ 
tions used in the program have been calculated by means 
of the same standard coalescence model as described 



FIG. 16: Atmospheric antiproton flux at mountain altitude 
(2770 m) measured by the BESS collaboration m (at the 
indicated geomagnetic latitude 0 m ), compared with the present 
calculations. The latter update the results presented in \6'tl . 
The grey histogram corresponds to the NAR process, the white 
histogram showing the full calculations. 


above, using the p production cross sections discussed in 
section ono In the propagation process, the absorption 
cross sections for p A collisions discussed in section !!!! A II 
above were used, as well as the NAR contributions for the 
p and d particles, discussed in section Hi! Bl 

Fig .IKilshows the p flux at mountain altitude (2770 m), 
i.e., deep inside the atmosphere, measured by the BESS 
collaboration El: compared to the present calculations. 
It can be observed that the ratio of the contribution of 
the tertiary p flux (grey histogram on the figure) to the 
parent secondary component, appears to be on the same 
scale as in the galactic flux, as expected. The overall flux 
appears to be somewhat («20%) smaller than the value 
reported in 0, due to the use of p A absorption (total 
reaction) cross section (see above) smaller in the present 
calculations than in the previous report. The larger rel¬ 
ative flux at small kinetic energies is due to the NAR 
process, taken into account in the present work while it 
was not in the previous one. The overall agreement with 
the data, although slightly less good than in 0, is still 
fair however. 

Fig. 1171 shows the downward flux of antimatter nuclei 
p , d , and t , at 400 km of altitude, induced by CR colli¬ 
sions with atmospheric nuclei, for a set of bins in latitude 
between equator and poles [72l ] . Similar distributions are 
obtained for the upstream flux, showing that the domi¬ 
nant part of this flux is made of trapped particles 0 . 
A striking feature of the results is the very large pro¬ 
portion of the flux due to the tertiary component (NAR 
process), while this fraction was much smaller both for 
the flux at ground level and for the galactic flux. This 
originates from the selectivity of the earth magnetic field 
which tends to trap more efficiently low energy particles, 
the NAR process decreasing the energy of the propagat¬ 
ing particle then trapped by the field and confined in the 
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belt. Note that in the polar region where the earth field 
is small, the whole flux is almost totally accounted for by 
the tertiary component. The same is true for the d flux 
(intermediate panels). The shape of the distributions 
of the secondary p and d fluxes are also found signifi¬ 
cantly different from their galactic counterpart (see |73l|l. 
In particular, the atmospheric secondary particle spectra 
appear to have their maximum at different energies than 
galactic particles. All these features are due to the dy¬ 
namics of the particles in the earth magnetic field which 
may induce considerable distortion of the primary parti¬ 
cle spectra [72• On Fig. (top panel), for the polar 
latitudes, the p flux compared with the TOA measure¬ 
ments by BESS shows that this atmospheric background 
is negligible at the considered satellite altitude for ener¬ 
gies T> 0.5 GeV. At lower energies however, the data 
and the atmospheric flux seem to converge towards close 
values. 

Fig. HH1 lower panel shows the calculated spectral dis¬ 
tributions for mass 3 and 4 antinuclei. The typical struc¬ 
ture has the same origin as for lighter antinuclei: the 
spectrum is limited by the rise of the particle production 
cross section driven by the production threshold (itself 
rising with the mass of the produced particle), on the 
lower side. The upper energy limit is set by the trap¬ 
ping conditions. For these fluxes the NAR contributions 
haven’t been included. Dynamically, the same consider¬ 
ations as given in section ITV Cl concerning the NAR pro¬ 
cess are valid as well for atmospheric particles. However, 
as seen above, the expected NAR effect is much larger 
for atmospheric than for galactic fluxes. Including this 
process would somewhat weaken the high energy compo¬ 
nent and strongly populate the low energies (empty on 
the figure). But the induced change (increase of the over¬ 
all flux), although significant at the considered order of 
magnitude, would not modify the conclusion, the overall 
flux being vanishingly small. 


VI. COUNTING STATISTICS FOR LIGHT 
ANTIMATTER NUCLEI AT SATELLITE 
ALTITUDE 


Particles p d 3 He 4 He 

Galactic rate 10 6 15 1CU 3 1CU' 

Atmospheric rate 5 10 5 3 10 -4 5 10~ 9 


TABLE III: Expected count rates for light antinuclei in the 
AMS experiment, for 3 years of effective counting MM , as 
described in the text. 


The forthcoming search experiments for primordial an¬ 
timatter will need some estimated values of the expected 
secondary galactic and atmospheric flux for the inves¬ 
tigated particles, to be confronted to the experimental 
results. The secondary fluxes of antimatter particles cal- 
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FIG. 17: Calculated fluxes of antimatter nuclei produced in 
the atmosphere, shown by groups of 4 panels, for p (top), 
d (middle), with the full (secondary+tertiary, white his¬ 
togram) and tertiary only (grey hi stog ram) fluxes shown. The 
recent experimental p data from \7C] for the galactic flux at 
TOA are also shown (full triangles). The lower panel shows 
the mass 3 and 4 antinuclei fluxes, t (dark grey), 3 He (white), 
and 4 He (light grey) (bottom). All calculations are for 400 km 
of altitude, and for the different bins in latitudes between equa¬ 
tor and poles indicated on each individual panel (0 m is the 
geomagnetic latitude angle in radian). 










































19 



FIG. 18: Top: Calculated atmospheric p flux at 38 km of alti¬ 
tude (top histogram) compared with the recent BESS98 mea¬ 
surements (which include both the galactic and the atmo¬ 
spheric components). The lower two histograms correspond 
to the calculated atmospheric d flux in the polar region, the 
gray histogram shows the tertiary component. The horizon¬ 
tal line labelled dsuSY corresponds to the limit for the d flux 
supersymmetric dark matter annihilation, from 0 / Bottom: 
Differential atmospheric p , d (x 5 1CP), 3 He (xl(f), fluxes 
in the polar region, at 38 km. For this latitude, the 3 He and 
t fluxes are identical (identical production cross sections). 


For the AMS experiment, a more accurate estimate 
of these numbers will be obtained with the processing 
of these fluxes in the AMS spectrometer simulation pro¬ 
gram, taking into account the identification capability of 
the instrument. The work is under way and the results 
will be reported later. 

VII. SUMMARY 

This work has provided new calculations of the light 
antinuclei p ,d , 3 He 4 He , fluxes produced by CR colli¬ 
sions on the ISM in the Galaxy, and the first calculations 
of the fluxes produced by CRs in the atmosphere for the 
same particles. For the d flux, a new production channel 
pp —► dX has been included in the calculations and was 
shown to contribute significantly to the total flux. The 
cross sections involved in the calculations were based to a 
large extent on existing accelerator data. The results are 
found significantly different from previous calculations. 
The calculated fluxes above 0.5-1 GeV/n are estimated 
to be reliable within a factor of 2 for the d flux, and within 
one order of magnitude for mass 3 antinuclei. The cal¬ 
culated flux in the low energy range, below 0.5 GeV/n, 
suffers larger uncertainties mainly because of the lack 
of experimental data. The calculated d flux appears to 
be marginally detectable by the AMS experiment mainly 
because it will have to be discriminated against a large 
background of ps and electrons. 

A further step will consist of using the calculated d flux 
as an input to the AMS02 spectrometer simulation to¬ 
gether with the other particles with the same electric 
charge (p , electrons) with their respective galactic fluxes, 
and discriminating the ds from the other particles by 
means of the spectrometer instruments, in order to eval¬ 
uate the capacity of the experiment to achieve the mea¬ 
surement of this flux. 

Acknowledgements: The authors are indebted to P. 
Salati and R. Taillet for several enlightening discussions 
on the subject. The diffusion model calculations are 
based on a code originally written by P. Salati. 


culated above have been energy integrated over the en¬ 
ergy range 0.1-100 GeV per nucleon and over the AMS 
orbit taking into account the effects of the geomagnetic 
cutoff and using an approximate acceptance of 0.5 m 2 sr 
for the AMS spectrometer. The results are given in ta¬ 
ble EH The resulting counting statistics to be expected 
is of the order of 10 6 , 20, 10 -3 , and 10” 7 particles per 
year, for p , d 3 He , and 4 He , respectively. Clearly, the 
rates for mass 3 and 4 antinuclei are negligibly small even 
at_the scale of very small statistics. Note however that a 
4 He candidate event would have the largest probability 
of being a misidentified background 3 He particle. 
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APPENDIX: INELASTICITY IN COMPOSITE 
PARTICLE SCATTERING 

Let us consider the system for incident deuterons 
rather than antideuterons, first. The result will be easily 
extended to d later. Qualitatively the effect of the small 
dissociation energy of the deuteron, or of any nucleus 
or other quantal system with small binding energy, is to 
make the collision highly diffractive because of the strong 
absorption induced by the large dissociation cross section 
of the deuteron. Because of this larger collision opacity, 
the inelastic collisions on nuclei will be more peripheral 
than they would be for an incident nucleon for example. 
For the same reason the angular distribution will be more 
forward peaked, partly because of the “fragility” of the 
projectile which limits the range of momentum transfer 
to low values, as explained below. But this does not im¬ 
ply that they would have a smaller cross section. In fact, 
it can be shown that to the contrary, the inelastic collision 
cross section induced by a composite nuclear projectile is 
expected to be larger than for incident nucleons. This is 
also supported by the existing data. 

In the case of dp collisions on a nucleon target however, 
the system is rather transparent and the diffractive effect 
can only be loosely referred to, the forward peaking of the 
cross section being induced mainly by the deuteron form 
factor. As discussed below, the latter is broader in the 
momentum space and the resulting angular distribution 
more forward peaked, which means that only low mo¬ 
mentum transfers are allowed in the scattering process. 

All the above arguments apply to the case of inci¬ 
dent antideuterons which collisions are made even more 
diffractive than for incident deuterons by the stronger 
absorption induced by the annihilation channels. The 
following arguments applies equally well to both cases. 

In the type of collision considered here, the incident 
d excites the ISM nucleon into a resonance-like state 
with mass Mj. The process is dominantly diffractive, 
characterized by small momentum transfer and periph¬ 
eral character. The hit nucleon can as well be bound 
inside a (ISM) nucleus. 

The transition amplitude F(t) (t 4-momentum trans¬ 
fer) for such a reaction can be most conveniently de¬ 
scribed for the present pedagogical purpose, in the single 
scattering limit of the Glauber multiple scattering ap¬ 
proach for the collisions of composite systems |X3, [7g, 
namely the impulse approximation, in terms of the form 
factors of the interacting systems G p (t ), G t (f) and of the 
interaction term v(t). The full Glauber amplitude can be 
written as : 

F{t) = j d 2 be< tf p #|T(&, a p , > (13) 

where 4>t{4>t) are projectile and target ground 


(excited) states wave functions respectively, b , s p , and 
St, the impact parameters variables, and T(b,s p ,st) = 
1 — e 1 ^ Xj the (Glauber) profile function of the collision 
system, with \j being the individual phase shift functions 
of the elementa ry c onstituents l74| , namely nucleons for a 
nuclear system |75j , constituent quarks for a subnucleon 
system (77ij, or both Q. This later relation is the foun¬ 
dation of the Glauber approximation and leads to the 
well known multiple scattering series of the scattering 
amplitude. 

A classical illustration of the correctness of the Glauber 
approach is precisely the application to the total deuteron 
cross section on proton, for which the above series for the 
case of elastic scattering using the optical theorem, leads 
to (7|: 

&tot(dp) ~ v(jpp) +<j{np ) — (2 nd_order.shadowingJ,erm) 

In the single scattering approximation, i.e., the Plane 
Wave Impulse Approximation (PWIA) (or Chou-Yang 
model for hadron-collisions), the elastic scattering am¬ 
plitude takes the very simple form: 

Fit) ~ A p A t G p {t)G t {t)v(t) (14) 

Where the numbers of constituents of the two systems A p 
and At originate from the normalization of the ground 
state wave function mm- From this relation, it can 
be shown straightforwardly by simple application of the 
optical theorem, that the total reaction cross section is 
proportional to the product of the number of constituents 
times the elementary total cross section between con¬ 
stituents a tot = A p At<TQ. This is the additive quark 
model for hadron-hadron collisions |80lj . 

The derivation for the inelastic Dp —> DX cross section 
goes through identical steps with the only difference that 
the proton final state is an excited state rather than the 
ground state, and the product of the initial and final pro¬ 
ton wave functions in the relation above which describes 
the ground state density in the case of elastic scattering, 
then becomes a transition density describing the excited 
nucleon state Sp(r) instead of the ground state density 

P O’)- 

A similar relation applies for inelastic scattering. In 
that case, one (or both) of the form factors involved de¬ 
scribe the excited state, namely: 

F[t) ~ A p A t G p {t)v(t)5 t (t) 

Therefore the expected inelastic Ap —> AX cross sec¬ 
tion, A being a composite system such as the deuteron d, 
is expected to be larger than its counter part in nucleon- 
nucleon scattering pp —> pX. 

This result is confirmed by the experimental data from 
f53| (see also |73|) where a larger cross section for dn —► 
d(pTT~) than for pn —> p(pTr~) was measured (see table 2 
and fig. 6 in [53|, see also 15711. 
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